What is in a pebble shape? 
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We propose to characterize the shapes of flat pebbles in terms of the statistical distribution 
of curvatures measured along the pebble contour. This is demonstrated for the erosion of clay 
pebbles in a controlled laboratory apparatus. Photographs at various stages of erosion are analyzed, 
and compared with two models. We find that the curvature distribution complements the usual 
measurement of aspect ratio, and connects naturally to erosion processes that are typically faster 
at protruding regions of high curvature. 

PACS numbers: 45.70.-n,83.80.Nb,91.60.-x,02.60.Jh 
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In soft matter physics it is well understood why cer- 
tain objects have naturally rounded shapes. For exam- 
ple, liquid droplets are round due to surface tension 
and bilayer vesicles are round due to bending rigidity [ 
In both cases the equilibrium shapes are determined by 
free-energy minimization, and sinusoidal perturbations 
decay exponentially with time according to a power of 
the wavevector. But how do pebbles in a river or on 
a beach come to have their familiar, smooth, rounded 
shapes (jjj? Are there similarly simple physical principles 
that govern the evolution of freshly-fractured polyhedral 
rocks into smooth rounded pebbles? Aristotle proposed 
that erosion is more rapid at regions farther from the 
center, where greater impulses can be more readily de- 
livered However this idea has not been quantified or 
tested, and we are aware of no other models for the origin 
of pebble shapes. 

A wealth of knowledge already exists on naturally- 
occurring pebbles JjL_ thanks to decades of field and labo- 
ratory study ■ There is a rich terminology 
for describing sizes (grain, pebble, cobble, boulder) and 
shapes (angular, rounded, elongated, platy). The sim- 
plest quantitative shape description is to form dimen- 
sionlcss indices from the lengths of long vs intermediate 
vs short axes. There is still debate as to how these axes 
are best defined, as well as to which dimensionless ratios 
are most useful [ill B EE Q EH • Another quantitative 
method used more for grains than for flat peb bles involves 
Fourier transform of the contour 0, IhH Il7l flj| . There 
the relative amplitudes of different harmonics give an in- 
dication of shape in terms of roughness at different length 
scales. However a Fourier description is not particularly 
natural because peaks erode faster than valleys, so that 
sinusoidal perturbations cannot erode without changing 
shape. To make progress in understanding the origin of 
pebble shapes, one needs a quantitative shape descrip- 
tion in terms of microscopic variables directly relevant to 
the erosion process. 

In this paper we propose to quantify the shape of fiat 
pebbles in terms of curvature, measured at each point 



along the entire two-dimensional contour. The curvature 
K equals the reciprocal radius of a circle that locally 
matches the contour, and can be deduced from the coor- 
dinates of the pebble boundary EH • This could be useful, 
since one can imagine faster erosion at protruding regions 
of higher positive curvature, and perhaps no erosion for 
scalloped regions of negative curvature. To illustrate the 
use of curvature for shape description, we perform a labo- 
ratory experiment in which clay pebbles are eroded from 
controlled initial shapes. We find that the curvature dis- 
tribution along the pebble contour becomes stationary, 
indicating a final fixed-point shape. Surprisingly, this 
stationary shape is not a perfect circle. We then intro- 
duce two models of pebble erosion, and test them using 
the curvature distribution as well as the aspect ratio. We 
find that the curvature distribution is a more incisive tool 
for discriminating between models. These results point 
the way for future studies, both of naturally-occurring 
pebbles and of lab-eroded pebbles of more direct geo- 
physical interest. 

For our laboratory experiments, the pebbles were pre- 
pared from ordinary white clay ( "chamotte" ) , molded 
into different polygonal shapes with uniform 0.5 cm 
thickness. Individual pebbles were eroded one at a time 
in a square pan, 30 x 30 cm 2 with 7 cm walls, spun at 
1 Hz around the central axis perpendicular to the bot- 
tom surface and oriented at 45° away from vertical. The 
pebble is first dragged upwards, at rest on the pan; near 
the top it begins to slide, and it accelerates until striking 
the wall; then it rolls along the side, comes to rest, and 
starts a new cycle. Erosion is mainly caused by collision 
with the walls. At five minute intervals, corresponding to 
roughly 300 collision events, the clay pebble is removed 
and photographed. As an example, photographs of one 
initially-square pebble with sides of 5 cm are shown in 
Fig-d digitized contours are shown in the inset of Fig.Qfe.- 
Evidently the erosion is fastest at the corners, which pro- 
trude and have high positive curvature. Once the corners 
have been removed, the pebble reaches a nearly round 
shape that progressively shrinks. 
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FIG. 1: Shape evolution of a 5 x 5 x 0.5 cm 3 square peb- 
ble eroded in a rotating basin, depicted on the lower right, 
with roughly 300 cutting events between each photograph. 
The time sequence is left to right, and top to bottom. The 
erosion is due to small random cutting events, a particularly 
large example of which is highlighted by arrow in the third 
photograph. 




FIG. 2: [Color online] (a) Digitized contours (inset) and cu- 
mulative curvature distribution for the initially-square clay 
pebble of Fig. Q Note that the red curves are indistinguish- 
able, indicating a stationary final shape. The final shape is 
not circular, as demonstrated both by the black-dotted circle 
in the inset and by the gradual rise of f(K) from zero to one. 
(b) Left axis: Cumulative curvature distribution for the fi- 
nal shapes attained by clay pebbles of varied initial shape, as 
labeled; the number of runs for each shape is given in paren- 
theses. Right axis: average curvature distribution for all final 
shapes, p — df/dK, and best fit to a Gaussian. 



The first step in analysis is to deduce the value of the 
curvature at each digitized point along the contour by 
linear-least-squares fits of radius vs angle to a third-order 
polynomial. Some care is needed in choice of fitting win- 
dow: if it is too small then the matrix inversion of the 
fit is singular j^jj, if it is too large then the fit deviates 
systematically from the contour points; results for all ac- 
ceptable windows are averaged together. As a check, the 
perimeter P is computed from straight-line segments be- 
tween adjacent points and is confirmed to equal 2-7r di- 
vided by the average curvature (K) 

As a shape-descriptor, we next construct the curva- 
ture distribution such that p(K)dK is the probability of 
finding a value of the curvature between K and K + dK. 
For reliability, it is preferable to work with the cumula- 
tive curvature distribution, f(K) = f Q p(K')dK', which 
may be computed as the fraction of the perimeter with 
curvature less than K [2(J. For a perfect circle, p{K) is 
a delta- function and f(K) is a step- function. The evolu- 
tion of f(K) for one square pebble is shown in the main 
plot of Fig. 12k, corresponding to the contours of the inset. 
Note that the results are plotted vs K/{K) in order to 
remove a trivial scale factor related to the ever-shrinking 
perimeter. Evidently, the curvature distribution starts 
broad but narrows down as erosion proceeds. By about 
the fifth (purple) or sixth (largest red) contour, the form 
of f(K) vs Kj (K) reaches a steady state that fluctuates 
around a well-defined average. If the initial shape is cir- 
cular, then by contrast the form of f(K) starts narrow 
but broadens as erosion proceeds. 

To test the possible influence of initial conditions on 
the final shape, we repeat the erosion experiment of 
Fig. [2k for clay pebbles molded into several different 
shapes. In all cases the form of f(K) vs K/(K) even- 
tually becomes time independent, with results given for 
different initial shapes on the left axis of Fig. [5Jd. To 
within fluctuations comparable with those seen for a sin- 
gle stationary pebble as it shrinks further, these are all 
identical. Thus the final shape is truly stationary, in- 
dependent of both size and the initial conditions, and 
may be quantified by averaging together the final cumu- 
lative curvature distributions for all pebbles; the results 
are shown by open circles in Fig. [2b- Note that its form 
is not set by uncertainty in curvature measurement, since 
the circular pebble started with a narrower distribution. 

The stationary curvature distribution, obtained by nu- 
merical differentiation of the average f(K), is shown on 
the right axis of Fig. 12b- Clearly it is not a delta-function, 
and hence the shape is not a perfect circle. The form of 
p(K) can be reasonably fit to a Gaussian, as shown, with 
a standard deviation of cr/(K) = 0.7. However, the ac- 
tual distribution is skewed in the tails and has a standard 
deviation of a/ (K) = 0.8. Though our laboratory erosion 
experiment may be artificial from a geophysical perspec- 
tive, it is interesting that a stationary shape exists and 
that it is not trivially circular. Clearly, f(K) was a useful 
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tool for this demonstration. 

In the remainder of the paper, we use the cumulative 
curvature distribution to evaluate the validity of two sim- 
ple models of the laboratory erosion experiment. We be- 
gin with a "polishing model" in which erosion is directly 
coupled to curvature. Namely, the normal velocity at 
each boundary point is taken in proportion to curvature, 
for regions of positive curvature, and is zero otherwise. 
(If all points are moved normally in proportion to cur- 
vature, irrespective of sign, then it has been proven that 
any initial contour will shrink to a point and approach a 
circle in this limit [HI-) The evolution predicted by the 
polishing model for a square pebble is shown in Fig. 
In the top plot, Fig.|2^, the contours and cumulative cur- 
vature distributions are plotted at the same sequence of 
perimeters, and with the same line-color codes, as the 
actual data in Fig.|3i. As expected, this model is highly 
efficient at polishing a rough pebble into a smooth shape. 
It becomes essentially circular, to the eye, soon after the 
entire original boundary has been eroded. The cumula- 
tive curvature distribution, f(K), appears to approach 
a step function, contrary to our laboratory experiments 
shown in Fig. [3 Thus the polishing model may be appro- 
priate for other natural or laboratory erosion processes, 
but not for ours. 

To track evolution toward a final shape more clearly, 
and to compare new vs old shape description methods, 
the aspect ratio of the pebbles and the width of the cur- 
vature distributions are shown in Figs. 0U>c. These are 
plotted vs perimeter divided by the initial perimeter, as 
a surrogate for the duration of erosion (except that time 
evolution is from right to left). Since major and minor or- 
thogonal axes are equal for a square shape, we work with 
the "caliper" aspect ratio, C/A, of largest to smallest sep- 
arations for parallel planes that confine the pebble; as the 
name suggests, this quantity may be readily measured in 
the field using calipers. For a dimensionless measure of 
the width of the curvature distribution, we take the stan- 
dard deviation, a, divided by the average curvature, (K). 
The results for all pebbles are plotted in Figs.|3>c, with 
four initially-square pebbles highlighted by larger sym- 
bols. Stationarity is reached when the square pebbles are 
eroded to about 3/4 of their initial perimeter, after which 
both C/A and o~/{K) data fluctuate significantly about 
an average value. These fluctuations represent real shape 
changes, not measurement uncertainty, and are larger for 
C/A than for a/ (K). Predictions by the polishing model 
for the evolution of the four square pebbles are also in- 
cluded for completeness. While the initial agreement is 
satisfactory, the polishing predictions for C/A and a/ (K) 
quickly fall below the data and, furthermore, do not ex- 
hibit fluctuations. 

To more successfully explain the laboratory erosion, 
we introduce a minimal one-parameter "cutting model". 
We imagine the effect of collision between pebble and 
wall is to create a straight fracture near the boundary 
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FIG. 3: [Color online] (a) Contours (inset) and cumulative 
curvature distribution predicted by the "polishing model" for 
a square pebble. The initial contour points, the sequence of 
perimeters displayed, and the line-color codes, are the same 
as in Fig. The open green circles represent the stationary 
f(K) averaged over all laboratory pebbles, (b) The caliper 
aspect ratio and (c) the standard deviation of the curvature 
distribution vs scaled perimeter; note that time evolution is 
from right to left. Symbol and line codes for (a) and (b) are 
as follows: Data for four different square pebbles are shown 
by large colored symbols; data for pebbles of other initial 
shape are shown as small gray pluses once stationarity has 
been reached; for the five pentagons P a was not measured, so 
data are shown arbitrarily at P/P = 0.1; aspect ratio and 
standard deviation for the final stationary shape of the labo- 
ratory experiments are shown by dashed lines. Predictions of 
the polishing model, applied to the four square pebbles, are 
shown by solid curves. 



whose size, on average, is set by the collision impulse; 
e.g. note the unusually large cut in the third photograph 
in Fig. ^ At each erosion step of the model, we thus pick 
a contour point at random and we choose a cut length 
from an exponential distribution whose average is some 
fraction a of the square-root of the pebble area. One can 
imagine other choices, but an exponential is simple and 
the impulse is set by the area rather than the perimeter. 

To match the cutting model with our data, we take 
a = 0.042 in order that the final, stationary shape have 
the same curvature distribution width as the laboratory 
pebble data. The resulting evolution of the contours 
and curvature distributions predicted for one square peb- 
ble are shown in Fig. 0Jl As observed in the lab, the 
corners erode faster than the flat sections and the final 
shape is nearly, but not perfectly, circular. The cutting 
model predictions are compared directly with the C /A 
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FIG. 4: [Color online] Predictions for evolution based on 
the "cutting model", and comparison with data. Displayed 
quantities are the same as for the "polishing model" in Fig. [3] 

and <j/{K) data in Figs. 0Jj-c. While the final average 
value of <j/{K) is correct by design, our model also cap- 
tures several other features of the data: the final average 
value of C/A, the magnitude of shape fluctuations in the 
stationary shape, and the rate of evolution from initial 
to final shapes. In spite of this success, the full form 
of f(K) in Fig. 0Ji is not completely correct in fine de- 
tail. Particularly, the model does not generate regions of 
negative curvature. No doubt this could be remedied by 
taking cuts from an appropriate distribution of shapes, 
not just straight. Such refinements could be tested by 
comparison with f(K) data. 

In conclusion, the distribution of curvature along the 
contour of a pebble is more incisive than the aspect ratio 
as a tool for both describing pebble shapes and for testing 
models. The characterization of such eroded forms is of 
obvious importance in the subfield of geology known as 
sedimentology Q, where pebbles are considered as wit- 
ness to the geological conditions under which they were 
formed. It would therefore be of great practical conse- 
quence if there existed a method to decipher the mes- 
sage imprinted in a pebble shape, a method that would 



distinguish beach vs glacier vs river erosion, that would 
tell how far has a pebble traveled down a stream and 
perhaps even for how long it has been subjected to ero- 
sion forces. Moreover, the ability to quantify and explain 
shapes resulting from abrasion, erosion or other wearing 
mechanisms would also be of significance in other fields, 
wherever nature or man transforms a solid object by ir- 
reversibly removing fractions of the mass in a sequence 
of elementary events. The curvature distribution could 
prove a helpful tool to study all such processes. 
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